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Abstract 

We introduce a new class of Fokker-Planck equations associated with an effec- 
tive generalized thermodynamical framework. These equations describe a gas of 
Langevin particles in interaction. The free energy can take various forms which 
can account for anomalous diffusion, quantum statistics, lattice models... When the 
potential of interaction is long-ranged, these equations display a rich structure as- 
sociated with canonical phase transitions and blow-up phenomena. In the limit of 
short-range interactions, they reduce to Cahn-Hilliard equations. 

Key words: Generalized Fokker-Planck equations, generalized thermodynamics, 
stochastic processes, long-range interactions 



1 Introduction 



The statistical mechanics of systems with long-range interactions is currently 
a topic of active research. So far, most works have focused on the case of 
Hamiltonian systems of particles such as the iV-stars problem, the iV-vortex 
problem or the HMF model [1]. For these systems the energy is conserved 
and the proper statistical description is the microcanonical ensemble. The 
statistical evolution of the particles is described by kinetic equations such as 
the Landau equation in astrophysics [2] or the kinetic equation derived in [3] for 
point vortices. Since statistical ensembles are not equivalent for systems with 
long-range interactions, it may be of interest, at a conceptual level, to compare 
this microcanonical evolution with a canonical one. In that respect, we can 
consider a system of Brownian particles in interaction. These particles are in 
contact with a thermal bath that imposes its temperature. Therefore, they 
have a rigorous canonical structure and their statistical evolution is described 
by non-local Fokker-Planck equations [4,5]. 
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In a different context, it has been shown that Tsallis generalized thermody- 
namics could be useful to interpret anomalous diffusion in complex systems 
and that the g-entropies are connected to nonlinear Fokker-Planck equations 
[6,7]. These equations can be derived from stochastic processes involving a 
multiplicative noise depending on the distribution function [8]. In fact, it is 
possible to construct a larger class of Fokker-Planck equations associated with 
a generalized thermodynamical framework encompassing Tsallis and Boltz- 
mann entropies [9-12,5]. These generalized Fokker-Planck equations are ob- 
tained from standard ones by assuming that the diffusion coefficient is an 
arbitrary function D(f) of the distribution function (the friction term can 
also depend on /). For these equations, there exists a Lyapunov functional, 
satisfying F < 0, that can be interpreted as a generalized free energy. In this 
context, Tsallis free energy corresponds to a power-law dependence of the dif- 
fusion coefficient D(f) ~ and the q parameter in Tsallis formalism is 
related to the exponent of anomalous diffusion. 

In this paper, following [11], we introduce a new class of Fokker-Planck equa- 
tions associated with a generalized free energy functional and an arbitrary po- 
tential of interaction. These equations connect the two topics discussed above: 
long-range interactions and generalized stochastic processes. Their mathemat- 
ical and physical richness is considerable, as they exhibit a wide variety of 
canonical phase transitions and blow-up phenomena [13-17]. In the limit of 
short-range interactions, they reduce to Cahn-Hilliard equations. 



2 Generalized Fokker-Planck equations in phase space 

2. 1 Langevin particles in an external potential 

We consider a system of non-interacting Langevin particles evolving in a fixed 
external potential $ ea;t (r). We assume that the motion of each particle is de- 
scribed by the generalized stochastic process 



where R(t) is a white noise satisfying (R(i)) = and (R a (t)Rb(t')) = 5 a b5(t — 
t') and C is an arbitrary convex function, i.e. C" > 0. Since the function in 
front of R(£) depends on (r,v), the last term in Eq. (1) can be interpreted 
as a multiplicative noise. We note, however, that it depends on (r, v) through 
the local density of particles /(r, v, t). This corresponds to a back-reaction 
from the macroscopic dynamics. When C(f) = ~zr[(f q — f), Eq. (1) reduces 



dr 
~dt 



dv 
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to the stochastic process studied by Borland [8] in connexion with Tsallis 
thermodynamics. 

Using standard methods, we can show that the stochastic process (1) leads to 
the generalized Fokker-Planck equation 

where we have introduced a generalized inverse temperature (3 = 1/T through 
the generalized Einstein relation £ = Df3. This equation can be put in the 
form 

where L is the advection operator and F is the generalized free energy 
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F[f] = E-TS = J f-d D vd D w + J p$ ext d D r + T j C(f)d D rd D v. (4) 



The first term in Eq. (4) is the kinetic energy, the second term is the potential 
energy and the third term is a generalized entropy S = — J C(f)d D rd D v . It 
is easy to show that the dissipation of free energy can be expressed as 

Vrd D v, (5) 



F = 



^~ fC"(f)^ + Pfv 



which is negative provided that D > 0. Therefore, the generalized Fokker- 
Planck equation (2) satisfies a canonical H-theorem F < 0. Hence, the free 
energy (4) is a Lyapunov functional. Finally, the stationary solutions of Eq. 
(2) are given by 

C*(/) = -0(y + (6) 



They extremize the free energy (4) at fixed mass and temperature. In addition, 
only minima of free energy are linearly stable via Eq. (2), see [11]. 



2.2 Langevin particles in interaction 



We now consider a system of Langevin particles in interaction described by 
the A^-body stochastic equations 
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Ri(0> (7) 



where /j = /(rj,Vj,i) and Rj(i) is a white noise acting independently on 
the particles. The particles interact via the potential U = J2i<j u ( r i — r j) 
where u(r,i — Yj) is an arbitrary binary potential depending only on the ab- 
solute distance between particles. When C(f) = —/In/, corresponding to 
the Boltzmann entropy, Eq. (7) describes a system of Brownian particles in 
interaction. 

Starting from the iV-body Fokker-Planck equation, using a Kramers-Moyal 
expansion and a mean-field approximation [4,5], we can derive from (7) the 
non-local generalized Kramers equation 



df df ^ df d . r , 

— + v • — - V$ • — = \D 

dt dr <9v <9v 



fC"(f)^ + Pfv 



(8) 



where $ is related to the density p = J fd D v by a relation of the form 

$(r,t) = j p{v',t)u{v-v')d D v'. (9) 



Equation (8) can be written as Eq. (3) where, now, the free energy is given by 



F[f] = J f^d D vd D v + l -J p$d D r + T J C(f)d D rd D v. 



(10) 



Furthermore, Eq. (5) remains valid so that F < 0. Finally, the stationary 
states of Eq. (8) are determined by the integro-differential equation 



C\f) = -/?{!£ + / /(r', v>(r - r')d D v'd D v} - a. 



(11) 



They extremize the free energy (10) at fixed mass and temperature. In addi- 
tion, only minima are linearly stable via Eq. (8), see [11]. 



2.3 The strong friction limit 

In the strong friction limit £ — > +oo, or equivalently for large times t ^> £ _1 , 
the distribution function is given by 



C'U) = -(3 



+ A(r,t) 



+ o(r x ), 



(12) 
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where A is related to the density p, using p = J fd D \ = p(A). Introducing the 
pressure p = jj f fv 2 d D v = p(A) and eliminating A, we find that the fluid is 
barotropic in the sense that p = p(p), where the equation of state is entirely 
specified by the function C(f). Furthermore, writing the hierarchy of moment 
equations [11] or using a formal Chapman-Enskog expansion [18], we can show 
that the evolution of the density is determined by the non-local generalized 
Smoluchowski equation 



dt 



|(Vp + pV$) 



(13) 



which can be written explicitly as 

p'(p)Vp + pV J u(r-r')p(r',t)d D r' 



dt U 



(14) 



This equation can be put in the form 



dt 



1 ^f SF \ 



(15) 



where F is the free energy 



F[p] = jpj P -^dp'd D v + 1J p(r, t)u(r - r')p(r', t)d D vd D v' . (16) 



It can be obtained from Eq. (10) by using the relation (12) valid in the strong 
friction limit £ — > +oo to express F[f] as a functional of p [11,18]. It is 
straightforward to show that 



F — — [ — (Vp + pV$) 2 d D r, 

J ip 



(17) 



which yields F < 0. Therefore, the free energy (16) is a Lyapunov functional for 
the generalized Smoluchowski equation (14). Finally, the stationary solutions 
of Eq. (14) are determined by the condition of hydrostatic equilibrium 



Vp + pV$ = 0, 



(18) 



which can be written explicitly as 



-Vp = -V j p(r>(r - r')d D r'. 



(19) 
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They extremize the free energy (16) at fixed mass. In addition, only minima 
are linearly stable via Eq. (13). 



3 Generalized Fokker-Planck equations in physical space 



3. 1 Langevin particles in an external potential 



We consider a system of non-interacting Langevin particles described by the 
generalized stochastic process 



j t =-iV<$> ext (v) + J2Dp 



C{ P ) 



R(t), 



(20) 



where £ is the drift coefficient and R(t) is a white noise. The evolution of the 
density is described by the generalized Fokker-Planck equation 



dp 
dt 



D 



P C"(p)Vp + (3pV$ ext 



(21) 



where we have introduced a generalized inverse temperature (5 — 1/T through 
the generalized Einstein relation £ = D(5. This equation can be put in the 
form 



dp 
dt 



5_F_ 
Jp~ 



(22) 



where F is the generalized free energy 

F[p] = E-TS = J p$ ext d D r + T J C(p)d D i 



(23) 



The rate of dissipation of free energy can be expressed as 
DT r i 2 



F = 



P 



pC"{p)Vp + (3 P V® ext 



al D r, 



(24) 



which is negative provided that D > 0. Therefore, the generalized Fokker- 
Planck equation (21) satisfies a sort of canonical H-theorem F < 0. Hence, 
the free energy (23) is a Lyapunov functional. Finally, the stationary solutions 
of Eq. (21) are given by 

C\p) = -P$ ext - a. (25) 
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They extremize the free energy (23) at fixed mass and temperature. In addi- 
tion, only minima of free energy are linearly stable via Eq. (21). 



3.2 Langevin particles in interaction 



We now consider a system of Langevin particles in interaction described by 
the iV-body stochastic equations 



dr / 
-^ = -SV i U(r 1 ,...,T N ) + J2Dp i 



Pi 



Ri(t), 



(26) 



where pi = p(rj, t) and R«(t) is a white noise. The potential of interaction U is 
defined as above. Starting from the iV-body Fokker-Planck equation, using a 
Kramers-Moyal expansion and a mean-field approximation [5], we can derive 
from Eq. (26) the non-local generalized Smoluchowski equation 



dp 
dt 



= V- \D 



pC"(p)Vp + /3pV$ 



(27) 



where <3> is related to the density p by a relation of the form (9). Equation 
(27) can be written as (22) where, now, the free energy is given by 



F = X - J p$d D r + T J C{p)d D r. 



(28) 



Furthermore, Eq. (24) remains valid so that F < 0. Finally, the stationary 
states of Eq. (27) are determined by the integro-differential equation 



C\p) = -pf P(r', t)u(r - r')d D r' - a. 



(29) 



They extremize the free energy (28) at fixed mass and temperature. In addi- 
tion, only minima are linearly stable via Eq. (27). 



3.3 Connexion with Cahn-Hilliard equations 



If we now consider the case of short-range interactions, it is possible to expand 
the potential 



$(r,t) = J u(t') P (t + r')d D r', 



(30) 
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in Taylor series for r' 



0. Introducing the notations 



a = 



u(\x.\)d D x. and b 




(31) 



we obtain to second order 



$(r,t) =ap(T,t) + -Ap(T,t). 



(32) 



In that limit, the free energy (28) takes the form 




+ V(p) }d D r, 



(33) 



where we have set V(p) = — (2T/b)C(p) — {a/b)p 2 . This is the usual expression 
of the Landau free energy. In general b is negative so we have to minimize this 
functional integral. For systems with short-range interactions, the conservative 
equation (27) becomes 



This is the Cahn-Hilliard equation which has been extensively studied in the 
theory of phase ordering kinetics. Its stationary solutions describe "domain 
walls" . We can view therefore Eq. (27) as a generalization of the Cahn-Hilliard 
equation to the case of systems with long-range interactions. 



4 Conclusion 

In this paper, we have introduced a wide class of generalized Fokker-Planck 
equations and we have shown that they were associated with a consistent ther- 
modynamical formalism. We have derived these equations from a particular 
stochastic process involving a special form of multiplicative noise, but they can 
also arise from other types of microscopic models. In our point of view, gen- 
eralized Fokker-Planck equations are "effective equations" trying to take into 
account "hidden constraints" in complex media that are difficult to formalize 
[12]. Since these equations decrease a functional similar to a generalized free 
energy, this leads naturally to a notion of "effective thermodynamics". This 
effective approach is not in contradiction with usual thermodynamics. This 
is just a heuristic approach trying to deal with complex situations in which 




(34) 
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standard thermodynamics is difficult to implement. We can either use ordi- 
nary thermodynamics (Boltzmann) and try to take into account microscopic 
constraints or keep only the most accessible constraints (mass, energy,...) and 
introduce an "effective entropy" . This heuristic approach may be useful in cer- 
tain occasions when an exact description of the system's dynamics is difficult 
to implement. Generalized thermodynamics is justified when all the accessible 
microstates are not equiprobable due to additional microscopic constraints. 
These constraints can be due to quantum mechanics (Pauli exclusion princi- 
ple) , to the existence of a lattice, to hard sphere effects, to a fractal structure 
of phase space, etc. We emphasize that our formalism is valid for more gen- 
eral functionals than the Tsallis entropy. Tsallis entropy forms, however, an 
important class of functionals associated with polytropic distributions and 
power-laws. These distributions generate a natural form of self-confinement 
that can be of interest in nonextensive systems. However, it was our interest 
here to show that a more general formalism could be developed consistently. 
A notion of "generalized thermodynamics" , with a different presentation and 
a different motivation, has been developed independently by Kaniadakis [9], 
Frank [10] and Naudts [19]. 

The non-local Fokker-Planck equations introduced in this paper can be of in- 
terest in different fields of physics, with various interpretations. In the usual 
case where they are associated with the Boltzmann entropy, they describe a gas 
of Brownian particles in interaction [4,5]. Therefore, they constitute the canon- 
ical counterpart of the microcanonical kinetic equations (Boltzmann, Landau, 
Lenard-Balescu,...) describing Hamiltonian systems of particles in interaction 
[2,3]. They are thus interesting at a conceptual level to test dynamically situ- 
ations of ensemble inequivalence for systems with long-range interactions. For 
example, the "isothermal collapse" of self-gravitating Brownian particles de- 
scribed by the Smoluchowski-Poisson system [13-16] can be compared with the 
"gravothermal catastrophe" of the Hamiltonian iV-stars system described by 
the Landau- Poisson system (or by the orbit averaged Fokker-Planck equation). 
Similarly, the Smoluchowski-Poisson system for Brownian point vortices can 
be contrasted with the kinetic equation obtained in [3] from the Hamiltonian 
dynamics. They exhibit a very different behaviour illustrating the inequiva- 
lence of statistical ensembles. 

At a practical level, the generalized Fokker-Planck equations presented in this 
paper and in [11] can be used as numerical algorithms to compute explicitly 
nonlinearly dynamically stable stationary solutions of the Vlasov and Euler 
equations. In that context, the generalized entropy S = — J C(f)d D rd D v is 
interpreted as a Casimir functional or a iJ-function [20]. Its maximization 
at fixed energy and mass (or circulation) determines a condition of dynami- 
cal stability, not a condition of thermodynamical stability. The fact that this 
variational principle looks similar to the variational principle that occurs in 
thermodynamics is essentially a thermodynamical analogy [11]. Since the gen- 
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eralized Fokker-Planck equations introduced in [11] increase the functional S 
while conserving E and M (or T), they can be used as physical algorithms to 
construct arbitrary nonlinearly dynamically stable stationary solutions of the 
Vlasov (or Euler) equation. Such stationary solutions can arise in practice as 
a result of a (possibly incomplete) violent relaxation, leading to a metaequi- 
librium state on the coarse-grained scale [21,11,12]. 

Finally, non-local drift- diffusion equations isomorphic to the Smoluchowski- 
Poisson system have a real interest in biology in relation with the concept of 
chemotaxis (see [17] for a description of the analogy between self-gravitating 
Brownian particles and bacterial populations). Now, it has been noted that the 
diffusion of bacteries is anomalous in general so that the diffusion coefficient 
which appears in the drift- diffusion equation depends on the density. This leads 
to a notion of "effective" generalized thermodynamics. This is simply because 
the true dynamics is so complex that it is replaced by a simple equation 
with an "ad hoc" diffusion coefficient whose specific form has to be deduced 
from experiments (this illustrates what we mean by "hidden constraints"). 
Accordingly, different forms of "generalized free energy" functionals can arise 
in biology depending on the context. 

In conclusion, the generalized Fokker-Planck equations presented in this pa- 
per and in [11] can be of considerable interest, in particular in the case of 
long-range interactions where they exhibit a rich variety of canonical phase 
transitions. A systematic study of this new class of equations has been initiated 
in [13-17]. 
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